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ABSTRACT 


There  are  several  issues  that  are  encountered  when  designing  a  Free  Electron  Laser  (FEL) 
oscillator  in  the  terahertz  (THz)  regime.  One  such  issue  is  that  there  are  few  inexpensive 
materials  that  are  suitable  for  forming  a  semi-transparent  out-coupling  mirror.  In  an 
attempt  to  mitigate  this  particular  issue,  the  focus  of  this  thesis  is  to  explore  various  out- 
coupling  techniques  for  a  notional  FEL  oscillator  in  the  THz  regime.  The  advantages  and 
disadvantages  of  several  out-coupling  methods  will  be  discussed,  including  hole  out- 
coupling,  and  double  mesh  out-coupling.  Simulation  results  will  be  presented,  showing 
the  effects  of  these  various  methods  on  FEL  extraction  and  optical  beam  quality. 
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I. 


INTRODUCTION 


A  Free  Electron  Laser  (FEL)  that  uses  electrons  that  are  not  bound  to  an  atom  to 
create  coherent  light.  This  type  of  laser  has  many  attributes  that  distinguish  it  from  the 
myriad  of  different  lasers  used  worldwide.  One  characteristic  is  that  FELs  are  tunable 
over  a  wide  range  of  the  electromagnetic  spectrum  and  can  be  designed  to  create  virtually 
any  wavelength  of  light,  from  microwave  to  X-ray.  This  is  because  FELs  do  not  depend 
on  a  conventional  gain  medium  to  produce  the  coherent  photons;  rather,  the  photons 
produced  are  a  result  of  the  electron  interaction  with  a  magnetic  and  optical  field.  The 
FEL  can  be  operated  at  much  higher  energies  than  other  conventional  lasers  since  there  is 
no  gain  medium  to  be  damaged. 

This  thesis  will  be  exploring  a  notional  FEL  operating  in  the  Terahertz  (THz) 
regime  of  the  electromagnetic  spectrum.  Operating  in  this  regime  does  present  some 
unique  issues.  The  first  issue  presented  is  that  of  clipping,  which  will  be  explained  in  a 
later  chapter.  Briefly,  clipping  is  an  issue  because  of  long  wavelength  associated  with 
THz.  Another  issue  is  that  many  materials  are  not  transparent  to  THz  light.  This  makes 
it  more  difficult  to  extract  usable  light  created  in  the  FEL  cavity.  This  issue  is  the 
motivation  behind  the  study  of  various  out-coupling  techniques  for  an  FEL. 

The  methods  of  out-coupling  to  be  explored  will  consist  of  hole  out-coupling  and 
mesh  out-coupling.  Hole  out-coupling  utilizes  a  highly  reflective  material  with  a  hole  in 
the  center  as  the  out-coupling  mirror  and  mesh  out-coupling  uses  a  fine  metallic  mesh 
grid  structure  in  place  of  the  mirror. 

The  notional  FEL  will  be  in  the  oscillator  configuration.  The  components  of  this 
configuration  will  be  described  later  in  this  chapter. 

A.  FEL  COMPONENTS 

A  FEL  in  the  oscillator  configuration  is  comprised  of  several  components.  The 
major  components  as  shown  in  Figure  1  include  the  injector,  the  accelerator,  the 
undulator,  the  optical  cavity,  the  energy  recovery  area,  and  the  beam  dump. 
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Figure  1.  FEL  in  the  oscillator  configuration.  Red  dots  are  electrons, 
purple  ellipses  are  light,  and  light  blue  ellipses 
are  out  coupled  light.  From  [1] 


1.  Injector 

The  injector  is  where  the  “free”  electrons  are  created  from  the  cathode.  The 
processes  by  which  electrons  are  emitted  from  the  cathode  surface  include  thermionic 
emission  and  the  photoelectric  effect.  Thermionic  emission  utilizes  a  heated  cathode  that 
emits  electrons  from  its  surface.  The  photoelectric  effect  uses  light  pulses  that  are 
incident  on  the  photocathode.  If  the  energy  of  the  light  is  greater  than  the  work  function 
of  the  cathode,  then  electron  pulses  will  be  emitted  [2],  The  emitted  electron  pulses  are 
then  accelerated  using  a  static  or  oscillating  electric  field  to  around  5  MeV. 

2.  Accelerator 

The  electron  pulse  next  enters  the  accelerator.  Many  accelerators  are 
superconducting  radio  frequency  linear  accelerators  (RF  LINAC)  that  have  niobium 
cavities  that  require  liquid  helium  for  cooling.  In  this  stage,  an  electron  pulse  travels 
through  each  cavity  in  phase  with  the  oscillating  RF  field  inside  the  cavity.  As  a  result, 
the  electrons  can  be  accelerated  to  a  highly  relativistic  energy,  100  MeV  for  typical 
applications.  The  FEL  can  have  electrons  that  have  energies  ranging  from  5  MeV  to  15 
GeV. 
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3. 


Undulator 


The  undulator  is  composed  of  a  series  of  permanent  magnets  that  are  arranged  in 
alternating  polarity  and  are  often  spaced  equally.  As  the  electrons  pass  through,  the 
Lorentz  Force  causes  them  to  “wiggle”  as  they  travel.  As  the  electrons  wiggle,  they 
create  a  light  pulse  that  has  a  wavelength  much  shorter  than  the  undulator  period  because 
of  relativistic  Doppler  shifts.  Only  a  small  percentage  of  the  energy  of  the  electron  pulse 
is  converted  into  light. 

4.  Optical  Cavity 

The  optical  cavity  is  the  region  of  the  FEL  where  coherent  light  is  amplified.  The 
components  of  this  region  are  only  used  in  the  oscillator  configuration.  The  oscillator 
configuration  uses  a  short  undulator  and  two  mirrors.  One  of  the  mirrors  is  highly 
reflective  and  reflects  the  created  light  pulse  back  through  the  undulator  for  stimulated 
emission  and  coherence  evolution.  The  out-coupling  mirror  is  typically  partially 
transparent  and  allows  for  a  small  percentage  of  the  light  to  pass  through. 

5.  Energy  Recovery  (The  Beam  Dump) 

Since  only  a  small  percentage  of  the  energy  of  the  electrons  is  converted  to  light, 
some  FELs  direct  the  electrons  to  pass  through  the  RF  accelerator  again  with  the  RF  field 
~  180-degrees  out  of  phase.  This  causes  the  electrons  to  lose  energy  to  the  RF  field,  and 
prevents  high  energy  electrons  from  entering  the  electron  beam  dump.  This  reduces  the 
amount  of  shielding  needed  for  any  applications.  This  type  of  configuration  is  call  an 
Energy  Recovery  LINAC  (ERL),  since  it  recovers  most  of  the  unused  energy  in  the 
electron  pulse,  thus  improving  the  overall  efficiency. 
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B.  FEL  THEORY 

1.  FEL  “Pendulum”  Equation  and  Electron  Dynamics 


The  dynamics  of  electron  motion  as  it  passes  through  the  undulator  are  governed 
by  the  following  equations: 


d{yp) 

dt 


e 

me 


(. E  +  fixB ) 


dy 

dt 


me 


(1.1) 


(1.2) 


r  = 


(1.3) 


where  E  and  B  are  the  electric  and  magnetic  fields,  respectively,  e  =  \e\  is  the  charge 
magnitude  of  an  electron,  m  is  the  mass  of  an  electron,  /3  =  v/c  is  the  dimensionless 

electron  velocity,  and  y  is  the  relativistic  Lorentz  factor  [1],  Equation  1.1  is  the  Lorentz 
force  equation  and  governs  the  motion  of  the  electron  as  it  encounters  FEL  fields. 
Equation  1.2  is  the  rate  of  energy  change  due  to  the  work  done  by  the  electric  field.  This 
equation  calculates  the  energy  given  to  the  optical  field  by  the  electron. 

As  the  electrons  travel  through  the  undulator,  they  interact  with  the  electric  and 
magnetic  fields  of  the  light,  E  =  £'(cos^,-sin^,0),  B  =  E (sin  y/,  cos  ^,0)  where  E  is  the 
field  amplitude  in  cgs  units,  y/ =  kz- cdt  +  tf),  with  (f)  the  optical  phase,  k  =  2njA  the 
wave  number,  A  the  wavelength,  and  co  the  optical  frequency.  The  electrons  also 
interact  with  the  magnetic  field  of  the  undulator,  B  =  5(cos(£0z),sin(&0z),o) ,  where  B  is 
the  undulator  field  strength,  and  A0  =  2 n/k^  is  the  undulator  period  [1].  The  electrons  at 
this  point  are  highly  relativistic,  so  fiz~  1 ,  so  when  the  fields  are  substituted  into 
Equation  1.1,  the  transverse  motion  in  the  undulator  is  approximately  determined  only  by, 
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(1.4) 


J3±  « - (cos(&0z),sin(k0z),o) 


r 

9 

where  K  -  eBiwsA()j^27imc2 )  is  the  dimensionless  undulator  parameter,  and  Brms  is  the 

rms  value  of  the  undulator  field.  Typically  K  «  1  for  most  FELs  [3].  If  Equation  1.4  is 
substituted  into  Equation  1.2,  the  electron’s  energy  evolution  yields 

eKE 

cos  + 

(1.5) 


y  =  —  =  — cos^-/?  sin^~|  =  ^^cos(<^+^) 
dt  me  L  y  J  ymc  v  ’ 


where  =(k  +  k^z- cotis  the  “electron  phase.”  The  electron  phase  follows  the 
microscopic  position  of  the  electron  as  it  evolves.  It  is  important  to  note  that  the  energy 
exchange  between  the  electron  and  the  optical  field  increases  when  cos(<^+  ^)  >  0  and 

decreases  when  cos(^  +  (f)  <  0 . 


The  “electron  phase  velocity”  is  defined  as  v  =  —£=L[[k  +  k0)/3,-k^  =  £, 

where  (...)  indicates  a  derivative  with  respect  to  dimensionless  time  t  =  ct/L  with  L  the 

length  of  the  undulator  so  that  T  varies  from  0  to  1  [1],  If  v=  0,  the  energy  exchange  is 
at  an  optimum  level  and  the  FEL  is  said  to  be  at  “resonance.”  In  this  condition,  the 
electron  phase  velocity  is  manipulated  to  show  that  the  resonant  laser  wavelength  is 


t  4>(i-A).,4>(i-*2) 
A  ~  27o2 


(1.6) 


where  yo  is  the  initial  electron  beam  Lorentz  factor  and  y0  »  1 . 


In  a  beam  of  electrons,  each  electron  will  have  a  unique  initial  phase 
do  )=(«=(k+Kh  as  they  enter  the  undulator.  In  a  typical  FEL,  ~106  electrons  are 

randomly  spread  in  each  section  of  the  electron  beam  ~/lin  length.  This  means  that 
roughly  half  of  the  electrons  will  lose  energy  to  the  laser  field  and  half  will  take  energy, 
as  described  in  Equation  1.5.  A  change  in  the  electrons  energy  will  also  create  a  change  in 

the  position,  ^and  the  electron  phase  velocity,  v.  The  change  in  the  phase  velocity  is 
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i>=  4;ry/y0  while  near  resonance  in  the  relativistic  limit  [1],  If  this  is  combined  with 
Equation  1.5,  then  the  FEL  “pendulum  equation”  is  derived 


v  =  4'  =  Hcos(4'+^)5  (L7) 

where  |a|  =  4 ftNeKLE / y^mc2  is  the  dimensionless  laser  field  amplitude[l].  From  this,  it 

is  shown  that  changes  in  the  design  of  the  FEL  will  affect  the  dimensionless  amplitude  of 
the  laser  field.  If  the  electron  beam  has  a  larger  Lorenz  factor,  the  strength  of  the  laser 
field  will  need  to  increase  in  order  to  create  the  same  dimensionless  amplitude. 

2.  FEL  Wave  Equation 

Having  solved  for  the  microscopic  motion  of  the  electrons  in  the  undulator,  the 
evolution  of  the  laser  field  in  the  undulator  must  be  derived.  Consider  the  full  vector 
wave  equation  with  the  Coloumb  gauge  [3] 


V- 


1  d2^ 
c2  dt2. 


(1.8) 


where  J±  =  — /?  £(3) 

[x-rt  (t)j  is  the  transverse  current  density,  r,  (?)  is  the  electron 

i 

position  at  time  t,  and  A^x,t )  is  the  optical  vector  potential.  The  electric  and  magnetic 
fields  of  the  laser  can  be  derived  using  the  vector  potential, 


E  =  -— _ 

c  dt  B  =  V  x  A 


(1.9) 


For  the  case  of  laser  light,  the  beam  is  partially  or  completely  coherent  and 
A  (x,t  j  has  a  slowly- varying  envelope  E(x,t )  in  the  direction  of  propagation  along  the  z 

axis  [3].  If  the  electric  field  is  assumed  to  be  complex  and  slow-varying  over  an  optical 
wavelength,  the  optical  vector  potential  is  assumed  to  have  the  form 


A 


^x,t)  =  ^^-sei[kz~at) 


(1-10) 
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where  E  =  \E\e"/>  is  the  complex  laser  electric  field,  and  £is  the  polarization  vector  of  the 


and  the  phase 


laser  field.  If  the  electric  field  amplitude 
be  slow-varying  along  the  z  axis  [3],  Equation  1.8  can  be  rewritten  as 

i(kz-6)t) 


x,t) 


are  assumed  to 


—*2  ^  .T  1 

r  d  15) 

Vi  +  2ik\ 

v  dz  c  dt) 

An  - 
E* — -Ji 
c 


an) 


Both  sides  of  the  equation  are  multiplied  by  ke  ,{t  <0^£  and  a  new  coordinate, 
u  =  z-ct,  is  introduced  in  order  to  invoke  the  method  of  characteristics!!].  With  that, 
Equation  I!  1  becomes  the  paraxial  wave  equation  with  a  source  current 


— -2  ^  i 

n  dY 

Vi  +  2  ik\ 

vc  dt) 

E  =  -—J1.£*e-i{kz-a,l) 


(1-12) 


Using  the  relationship  between  the  transverse  motion  of  an  electron  and  the 
transverse  current  density  defined  previously,  Equation  I!  2  becomes 


—2  2  ik  d 

v±+ — — 

c  dt 


E  =  - AmeKkp{x,t )  <  e  ^  jy 


>(P) 


(1-13) 


where  pyx,tj  is  the  local  electron  density  in  a  small  volume  element  dV located  at  yx,tj, 
and  <. .  .>  represents  an  average  over  sample  electrons  in  that  volume  element. 

Dimensionless  time,  t,  is  reintroduced  and  the  wave  equation  is  multiplied  by  the 
factor  -AnNeKL2 / ylmc2k  so  that  it  becomes 


iL  -2  d 

- V_|_  H - 

2k  dr 


(1.14) 


where  a  =  \a\e"f'  is  the  complex  laser  field,  \ci\  =  AjiNeKLEj ypnc1  is  again  the 

dimensionless  laser  field  amplitude,  j  =  $7r2Ne2K2L2 pj y30mc2  is  the  “dimensionless  FEL 
density,”  and  it  is  assumed  that  y  ~y0  for  all  the  electrons  in  the  beam[3]. 
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The  wave  equation  is  made  completely  dimensionless  by  defining  the 
dimensionless  transverse  coordinates  x~x{k/2L)  ,  anc[  >'  The 

equation  is  now  written  as 


wave 


l  =2 


--V1  +  — 
4  dr 


(1.15) 


where  the  coordinates  (x,  y,  z,  r)  are  now  all  dimensionless.  When  the  laser  beam  and  the 
electron  beam  overlap  and  diffraction  is  small  over  the  undulator  length,  the  FEL  wave 
equation  can  be  written  in  its  simplest  form 


a  =  -j  <e  ^  > 


(L16) 
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II.  THZ  ISSUES 


A.  MATERIAL  CONSIDERATIONS 

For  the  purpose  of  this  thesis  THz  radiation  is  defined  to  be  in  the  range  of  100 
microns  to  1mm.  There  are  several  issues  that  arise  when  out-coupling  THz  radiation  in 
a  FEL  oscillator  in  this  range.  Typically,  FELs  use  a  partially  transparent  mirror  to  out- 
couple  the  light  from  inside  the  optical  cavity.  These  mirrors  consist  of  a  substrate  that  is 
then  layered  with  a  coating  that  is  reflective  to  the  wavelength  of  light  being  used.  As 
such,  there  are  several  material  properties  of  the  mirror  substrate  that  must  be  accounted 
for  in  order  to  produce  usable  THz  light  for  the  FEL. 

1.  Transmission  of  Radiation 

In  order  for  a  FEL  to  operate,  a  fraction  of  the  light  produced  in  the  optical  cavity 
must  be  transmitted  through  the  out-coupling  mirror  of  the  resonator.  THz  radiation  is 
absorbed  by  many  materials,  which  often  leads  to  the  need  for  expensive  materials  to  be 
used. 

The  fraction  of  incident  radiation  absorbed  can  be  calculated  by: 

A  =  \-{\- R)e~ax 

where  A  is  the  fraction  of  radiation  absorbed,  R  is  the  reflection  coefficient,  a  is  the 
absorption  coefficient,  and  x  is  the  thickness  of  the  material  [4],  This  shows  that  with  an 
increase  of  the  thickness  of  the  material  or  with  an  increase  in  the  absorption  coefficient, 
the  fraction  of  light  absorbed  increases.  This  absorbed  energy  manifest  itself  as  heat  and 
can  result  in  other  effects  on  the  mirror. 


2.  Thermal  Expansion 

When  a  majority  of  materials  are  heated  by  a  source,  there  is  a  tendency  for  the 
materials  to  expand.  This  thermal  expansion  is  found  experimentally  and  calculated  by: 


aT 


L 


dL\ 
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where  aj  and  f>  are  the  coefficients  of  linear  and  volumetric  expansion,  respectively,  L 
and  V  are  the  materials  length  and  volume,  respectively,  and  T  is  temperature  [5].  From 
these  two  equations,  it  can  be  seen  that  in  order  to  have  a  material  that  exhibits  minimal 
change  in  dimension,  the  material  in  question  must  have  a  small  coefficient  of  expansion. 
A  minimal  change  in  dimension  is  critical  in  order  to  prevent  optical  distortions,  which 
can  lead  to  the  laser  either  not  operating  at  the  optimal  levels  or  not  operating  at  all. 


3.  Thermal  Conductivity 

Thermal  conductivity  describes  the  ability  of  materials  to  conduct  and  transfer 
heat.  Thermal  conductivity  is  found  experimentally  and  calculated  by 

k=  L  fdQ ' 

dt  j 

where  k  is  the  thermal  conductivity  in  watts  per  meter-Kelvin,  L  is  the  lateral  distance 

between  the  temperature  difference  AT,  and  is  the  rate  of  heat  flow  into  the 

material.  [6]  The  thermal  conductivity  is  shown  to  be  proportional  to  the  heat  flow.  As 
such,  a  material  with  a  large  k  value  will  be  able  to  remove  heat  better  than  one  with 
small  values. 

Because  the  out-coupling  mirror  absorbs  large  amounts  of  energy  from  the  FEL,  it 
is  important  that  the  material  that  comprises  the  lens  be  able  to  transfer  the  heat  from  the 
focal  point  of  the  beam.  If  this  is  not  achieved  the  beam  will  distort  and/or  potentially 
destroy  the  mirror. 


4.  Insufficiencies 

There  have  been  some  studies  and  experiments  conducted  to  identify  the  optimal 
material  for  use  with  THz  radiation.  Table  1  provides  a  summary  of  selected  materials 
with  experimental  data  for  the  previously  discussed  properties.  The  color  scheme  of  the 
table  is  used  to  highlight  material  property  values  in  a  range  that  would  be  suitable  for 
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use  in  a  FEL  operating  in  the  THz  range.  As  shown,  only  one  of  the  five  materials, 
Chemical  Vapor  Deposition  (CVD)  Diamond,  appears  to  be  suitable  for  use  as  a  mirror. 
However,  even  CVD  Diamond  has  its  drawbacks.  Because  CVD  Diamond  is  very 
expensive,  this  would  lead  to  a  very  expensive  mirror  regardless  of  the  dimension. 


Materials 


Thermal 
Conductivity 
@20  °C  (W/M  °C) 


Thermal 
Expansion 
Coefficient 
(xlO-6  /  °C) 


Transmission 

Range 

(pm) 


HRFZ  Si 

1.5 

2.6 

50-300 

Quartz 

1.3 

0.59 

Sapphire 

42 

5 

Zinc  Selenide 

18 

7.57 

10-20 

CVD  Diamond 

3300 

1.2 

10  -  300 

Table  1.  Summary  of  selected  materials  and  respective  properties.  From  [7] 


B.  CLIPPING 

FEL  Beam  Clipping  is  the  process  by  which  an  optical  beam  is  negated  when  it 
comes  in  contact  with  an  absorbing  material.  Because  of  this  interaction,  a  portion  of  the 
beam’s  energy  is  lost  which  results  in  a  decrease  in  power.  This  phenomenon  is  most 
prevalent  in  the  undulator  section  of  the  FEL  because  the  beam  line  is  constricted. 
Figure  2  illustrates  how  a  clipping  interaction  in  the  undulator  would  occur. 


Figure  2.  Clipping  in  the  undulator  as  a  function  of  diffraction. 
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A  contributor  to  the  degree  of  severity  of  the  clipping  that  occurs  is  beam 
diffraction.  Beam  diffraction  can  be  described  by 


w(z)  =  w0 


(II.  1) 


where  w(z)  is  the  radius  of  the  beam  at  distance  z,  wo  is  the  optical  mode  radius  and  zR  is 
the  Rayleigh  length  [1],  From  Equation  II.  1,  as  the  distance  from  z=0  increases  the  radius 
of  the  beam  also  increases.  This  effect  is  further  exaggerated  because  the  Rayleigh 
length,  defined  by 


n  nw; 


ZR=- 


0 


/l 


(II.2) 


where  A  is  the  optical  wavelength,  and  n  ~1  is  the  index  of  refraction  decreases.  As  the 

wavelength  increases  the  Rayleigh  length  increases  which  in  turn  increases  the  beam 
radius  at  distance  z.  This  causes  the  beam  to  expand  rapidly  and  as  such  increases  the 
amount  of  the  beam  that  potentially  will  be  clipped. 

The  potential  for  beam  clipping  is  very  prevalent  with  THz  radiation  and  must  be 
taken  into  account  in  accordance  with  previous  studies  [7]  .As  such;  the  simulations  that 
will  be  conduct  will  use  a  constant  clipping  factor  of  4%  per  pass  through  the  undulator. 
A  separate  study  of  how  clipping  affects  the  extraction  of  energy  of  the  FEL  operating  in 
the  THz  regime  will  also  be  conducted  in  a  later  chapter. 
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III.  GAUSSIAN  MODES 


A.  GAUSSIAN  BEAMS 

There  are  several  ways  to  model  the  light  wave  produced  by  an  FEL  oscillator. 
The  most  convenient  approach  is  to  use  Gaussian  modes.  A  Gaussian  beam  in  the 
fundamental  transverse  mode  describes  a  light  wave  where  the  power  and  intensity  of  the 
light  are  at  a  maximum  at  the  center  of  the  beam  and  decrease  as  the  distance  from  the 
center  increases.  This  approximation  is  a  solution  to  the  paraxial  wave  equation. 

1.  Short  Derivation  of  Paraxial  Wave  Equation 

Consider  the  wave  equation  found  by  solving  Maxwell’s  equations  in  free  space: 

v2£=«02%2;  (mi) 

where  the  spatial  derivatives  are  taken  to  be  in  rectangular  coordinates.  If  it  is  assumed 

that  E  is  polarized  in  one  direction,  so  that  the  scalar  form  of  the  wave  equation  can  be 
invoked 

=  (1,1.2) 

where  l/ c 2  =  ju()s() .  It  can  be  seen  that  E  has  a  spatial  and  temporal  dependence.  As 

such,  the  electric  field  can  be  defined  as  the  product  of  both  spatial  and  temporal  factor  so 
that 

E(x,  y,  z,  t)  =  E{r ,  t)  =  E(r)eiat ;  (III.3) 

where  co  =  (2/re)  j  X  =  kc .  When  substituted  into  III.2,  the  wave  equation  becomes 

[v2  +  k2)E  =  0;  (III. 4) 

where  k 2  =  of  j  c1  ,and  the  explicit  time  dependence  has  been  removed  from  the  wave 
equation  [8]. 
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To  arrive  at  the  paraxial  wave  equation  from  Equation  III.4,  it  is  assumed  that  the 
transverse  beam  profile  changes  slowly  on  the  scale  of  an  optical  wavelength.  As  such, 
we  assume  the  form 

E(x,y,z)  =  u(x,y,z)e~,h-,  (III.5) 


where  u  describes  the  transverse  profile  of  the  traveling  beam  [8].  In  other  words,  u 

changes  slowly  in  z  relative  to  the  factor,  e~,kz .  When  Equation  III.  5  is  substituted  into 
III.4  the  wave  equation  becomes 


a2  ~  d2  ~  d2  ~  d  ~ 

— -u-\ - -u-\ - -u-2ik — u  =  0. 

dx2  dy2  dz2  dz 


(III.6) 


In  III. 6  the  second  derivative  of  the  z  component  of  u  can  be  assumed  to  be  negligible  by 
invoking  the  paraxial  approximation,  32«/<9z2|  <§:  \2kdu/dz\  [8],  From  this,  Equation  III. 6 
can  be  rewritten  as  the  paraxial  wave  equation: 


f 


d 


V  -2 ik—, 
1  dz) 


u  =  0; 


(in.  7) 


where  V2  =  d2 /dx2  +  d2/dy2 .  This  approximation  remains  valid  for  traveling  beams  that 

propagate  at  an  angle  less  than  30  degrees  with  the  z-axis.  Due  to  this,  the  paraxial 
approximation  holds  true  for  most  FELs  in  the  oscillator  configuration. 

B.  HERMITE-GAUSSIAN  BEAMS 

A  solution  to  the  paraxial  wave  equation  must  be  found  in  order  to  properly  model 
a  propagating  beam.  One  such  solution  is  the  Hermite-Gaussian  solution, 


Elm(x,y,z)  =  E0-^Hl 
w(z) 


(  \ 
w{z)J 


H„ 


f  \ 

V2  y 


v 


w(z) 


) 


xexp 


-ik  -  +~^  -  X  — ikz  +  i(l  +  m  +  l)^(z) 


27?  (z)  w2  (z) 


(III.8) 
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This  solution  is  very  intricate  yet  it  is  still  in  the  form  of  E\x,y,z )  =  u[x,y,z^e~lkz 
[8],  The  waist  size  term,  w(z),  was  initially  defined  by  Equation  II.l.  R(z)  is  the  radius  of 
curvature  of  the  propagating  beam  and  it  is  given  by  R[z}  =  z[\  +  [z2R/  z2^.  From  this 

definition,  it  can  be  seen  that  the  radius  of  curvature  changes  depending  on  whether  the 
beam  is  propagating  in  the  positive  or  negative  direction.  rj(z)  represents  a  phase  shift 

defined  by  ^^zj  =  arctan^z/zflj.  The  /  and  m  labels  are  integer  indices  of  the  Hermite 

polynomials,  Hi  and  Hm.  These  indices  ultimately  affect  the  transverse  intensity  of  the 
beam  and  are  referred  to  as  modes.  The  1=0,  m=0  mode  is  called  the  “fundamental” 
mode  because  the  subsequent  Hermite  polynomials  are  both  more  structured.  Figure  3 
displays  several  of  the  lower  order  spatial  modes  that  the  beam  can  assume. 


Figure  3.  Transverse  Electromagnetic  Modes.  From  [9] 

For  the  purpose  of  the  FEL,  the  fundamental  mode  is  often  preferred  because  the 
beam  can  be  later  focused  into  a  tighter  spot  for  a  given  set  of  focusing  elements.  For  the 
purposes  of  directed  energy,  this  means  more  power  can  be  delivered  to  a  smaller  spot  on 
a  distant  target. 
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IV.  HOLE  OUT-COUPLING 


With  hole  out-coupling,  the  light  created  in  the  FEL  optical  cavity  is  transmitted 
through  a  hole  in  a  highly-reflecting  mirror  instead  of  being  transmitted  through  a 
partially  transparent  mirror.  This  technique  has  been  shown  to  work  for  other  types  of 
lasers  operating  at  various  wavelengths  [10]. 

A.  DESIGN  PARAMETERS 

For  the  hole  out-coupling  simulations,  Table  2  lists  the  design  parameters  that 
were  held  constant: 


Eb  ~  5  MeV 

f  ~  50  MHz 

N=  16 

Xq  ~  2  cm 

X  ~  100  pm 

K  ~  0.35 

Table  2.  FEL  Parameters 


where  Eb  is  the  beam  kinetic  energy,  f  is  the  pulse  repetition  frequency,  N  is  the  number 
of  undulator  periods,  Xq  is  the  undulator  period,  X  is  the  optical  wavelength,  and  K  is  the 
dimensionless  undulator  parameter. [3] 

B.  SIMULATION  RESULTS 

Simulations  were  conducted  using  the  NPS  3D  oscillator  FEL  code  and  were 
categorized  into  two  groups  based  on  the  electron  bunch  charge,  lOOpC  and  the  200pC. 

1.  lOOpC  Bunch  Charge  Simulations 

Simulations  were  conducted  for  the  lOOpC  bunch  charge  varying  the  radius  of  the 
hole,  the  quality  factor  Q,  and  the  clipping  gap. 

a.  Varying  Hole  Radius  with  a  high  0  value 

In  the  first  set  of  simulations  conducted  for  the  lOOpC  bunch  charge,  the 
hole  radius  was  varied  over  a  range  of  values.  The  quality  factor,  Q,  is  defined  to  be 
inversely  proportional  to  the  amount  of  loss  in  a  cavity  without  a  hole.  For  this  set  of 
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simulations,  Q  was  set  at  a  very  large  value  (~1010)  in  order  to  approximate  a  highly- 
reflecting  mirror  surface.  This  means  that  essentially  no  light  is  absorbed  by  the  mirrors, 
thus  the  only  way  for  the  light  to  escape  is  through  the  hole.  The  hole  radius  range  was 
from  0.05  cm  to  0.40  cm.  The  simulation  results  are  plotted  in  Figure  4.  The  shape  of 
the  plot  shows  that  there  is  a  peak  of  -2.8%  extraction  at  a  hole  radius  of  0.125  cm.  This 
is  due  to  how  the  electrons  are  evolving  to  produce  the  light,  as  will  be  seen  in 
subsequent  figures. 


Extraction  for  a  Highly  Reflective  Surface  with  a  Hole 


Figure  4.  Simulations  results  for  lOOpC  bunch  charge  with  Q=  1010 

Figure  5  (a)  displays  the  optical  field  evolution  in  the  undulator  for  the 
hole  radius  of  0.05  cm.  The  blue/cyan  area  represents  the  optical  field,  the  red  dots  are 
the  electrons,  and  the  yellow  curves  represent  the  mode  of  the  beam  at  the  ends  of  the 
undulator.  In  this  case  the  peak  dimensionless  optical  amplitude  is  high,  |a| « 52 . 

indicating  FEL  saturation  in  strong  optical  fields..  With  a  small  hole  only  a  small 
fraction  of  the  light  produced  in  the  optical  cavity  is  out-coupled.  The  optical  field  grows 
to  saturation,  while  the  electrons  also  evolve  described  by  the  “pendulum”  equation  and 
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the  FEL  wave  equation.  In  Figure  5  (b),  the  electrons  appear  to  overbunch  as  they  travel 
farther  in  their  phase  space  orbits,  absorbing  energy  from  the  optical  field  and  hence 
reducing  the  extraction. 


Figure  5.  Simulation  output  for  hole  radius  of  0.05  cm.  (a)  The  optical  field  evolution  in 
the  undulator.  (b)  The  electron  evolution  phase  space  plot  at  the  end  of  the  undulator. 


The  optimal  hole  radius  of  0.125  cm  gives  the  largest  amount  of  optical 
power  out  of  the  FEL  mirror  hole.  The  peak  optical  field  is  smaller  ( |a|  =  27 )  than  that  of 

the  smaller  hole  radius,  since  a  greater  fraction  of  the  light  is  out-coupled  per  pass.  Figure 
6  (b)  shows  the  electrons  do  not  overbunch  as  severely,  and  hence  more  energy  is 
converted  to  light  and  out-coupled  through  the  mirror. 


0  t  1 


(a)  (b) 

Figure  6.  Simulation  output  for  hole  radius  of  0. 125  cm. 

If  the  hole  radius  becomes  too  large,  too  much  of  the  light  escapes  from 
the  cavity  which  causes  a  drop  in  the  dimensionless  field  amplitude  (|a|  =  2.4)  as  seen  in 
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Figure  7  (a).  Due  to  this,  the  electrons  do  not  evolve  much  in  phase  space,  as  seen  in 
Figure  7  (b),  and  very  little  energy  is  extracted  from  them. 


(a)  (b) 

Figure  7.  Simulation  output  for  hole  radius  of  0.25  cm. 


The  use  of  hole  out-coupling  did  have  a  negative  effect  on  the  transverse 
mode  shape.  As  seen  in  Figure  8,  as  the  size  of  the  hole  increased,  small  perturbations  in 
the  mode  shape  at  the  out-coupling  mirror  started  to  develop.  The  beam  is  losing  energy 
to  the  formation  of  higher-order  modes,  which  is  not  a  favorable  condition  for  the  FEL 
interaction  and  applications. 


a  (x,  y,  5 . 5)  6  .  l|a  (x,  y,  5 . 5)  3 . 6]  a  (x,  y ,  5 . 5) 0 . 3 

^  ^ 


(a)  (b)  (c) 

Figure  8.  Transverse  mode  shape  at  the  out-coupling  mirror  for  a  hole  radius  of:(a)  0.05  cm, 

(b)  0.125  cm,  and  (c)  0.25  cm. 

b.  Varying  Hole  Radius  with  a  Q  value  of  50 

For  this  set  of  simulations,  the  quality  factor  was  Q= 50,  representing  a 
loss  of  1%  of  the  energy  due  to  absorption  at  each  mirror  surface.  The  hole  radius  was 
again  varied  from  0.05  cm  to  0.40  cm.  The  range  was  selected  to  provide  a  direct 
comparison  between  the  highly-reflecting  mirror  surface  and  the  more  realistic  mirror 
surface. 
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Figure  9  plots  the  overall  results  of  the  simulations  showing  a  trend 
towards  smaller  hole  radii.  A  smaller  hole  radius  is  expected  to  be  necessary,  since  now 
that  some  of  the  energy  is  absorbed  by  the  mirror  surfaces,  the  out-coupling  hole  should 
be  smaller  to  compensate  for  the  absorption.  Due  to  limitations  on  the  numerical  grid,  the 
optimal  hole  radius  was  unable  to  be  found,  which  appears  to  be  less  than  0.05  cm  for 
this  case. 


Extraction  for  a  Partially  Absorbant  Surface  with  a  Hole 


Figure  9.  Simulations  results  for  lOOpC  bunch  charge  with  Q=  50 

c.  Varying  the  clipping  gap  for  the  theoretical  simulations 

Since  the  longer  wavelengths  expand  more  rapidly  than  shorter 
wavelengths,  absorbers  were  simulated  at  the  ends  of  the  undulator  to  clip  a  portion  of  the 
field  that  might  scrape  in  the  ends  of  the  undulator.  The  gap  between  the  absorbers  was 
varied  from  0.3  cm  to  3.0  cm  while  the  hole  size  was  kept  constant  at  0.125  cm.  For  the 
results  in  the  previous  section,  the  absorbing  gap  was  kept  constant  at  1 .4  cm. 

Figure  10  shows  the  results  of  the  simulations,  and  finds  that  there  is  an 
optimal  clipping  gap  of  1.5  cm  for  this  set  of  parameters.  As  the  gap  size  increases  past 
the  optimal  there  is  little  change  in  the  extraction  levels.  This  intuitively  makes  sense 
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because  if  the  size  of  the  gap  keeps  increasing,  eventually  there  is  no  clipping.  Past  that 
threshold,  there  should  be  no  change  to  the  extraction  since  no  energy  is  being  lost  to  the 
absorbers. 


Gap  Varitions  for  Bunch  Charge  of  100  pC 


Figure  10.  Simulation  results  for  Gap  variations. 


Figure  11  shows  the  simulation  output  for  the  optimal  gap  of  1.5  cm.  It 
can  be  seen  in  that  the  mode  shape  is  starting  to  become  slightly  less  Gaussian  in  (a)  and 
(c).  This  means  that  slightly  more  energy  is  going  into  higher-order  modes,  and  shows 
that  clipping  can  be  used  to  force  the  optical  mode  into  the  fundamental. 


Figure  1 1 .  Simulation  output  for  gap  of  1 .5  cm. 
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2.  200  pC  Bunch  Charge  Simulations 

Simulations  were  also  conducted  for  the  200  pC  bunch  charge,  again  varying  the 
same  parameters:  the  quality  factor  Q,  the  radius  of  the  hole,  and  the  clipping  gap. 

a.  Varying  Hole  Radius  with  a  high  Q  value 

In  the  same  way  as  with  the  100  pC  simulations,  the  set  of  simulations 
were  conducted  over  a  range  of  hole  radii  with  Q  ~1010.  The  same  hole  radius  range  of 
0.05  cm  to  0.40  cm  was  utilized. 

The  results  of  the  simulations  are  shown  in  Figure  12.  The  results  show 
an  optimal  value  for  the  hole  radius  of  0.15  cm  with  a  corresponding  extraction  of  2.9%. 
The  optimal  hole  size  is  slightly  larger  than  the  100  pC  case  since  there  is  more  gain  per 
pass  with  the  200  pC  bunch,  which  allows  more  losses  in  steady-state. 

Extraction  for  a  Highly  Reflective  Surface  with  a  Hole 


Figure  12.  Simulations  results  for  200  pC  bunch  charge  with  Q=  1010 

Figure  13  displays  the  output  results  for  the  optimal  hole  radius  of  0.15  cm 
showing  a  well-bunched  distribution  in  phase  space  in  Figure  13  (b)  and  slight 
perturbations  in  the  mode  shape  due  to  the  hole  in  Figure  13  (c). 
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b.  Varying  Hole  Radius  with  a  Q  value  of  50. 

For  this  set  of  simulations,  the  quality  factor  was  Q  =  50.  As  previously 
mentioned,  this  represented  a  loss  due  to  absorption  at  the  mirror  surfaces.  The  hole 
radius  range  again  was  from  0.05  cm  to  0.40  cm. 

Figure  14  illustrates  the  overall  results  of  this  set  of  simulations  directly 
compared  to  the  100  pC  case.  It  shows  that  there  was  an  increase  in  the  optimal  hole 
radius  compare  to  the  smaller  value  as  seen  in  the  previous  case.  For  this  set  of 
simulations,  the  optimal  hole  radius  is  0.10  cm  and  produced  an  extraction  of  2.9%.  The 
optimal  hole  radius  is  larger  than  the  100  pC  case,  since  the  higher  gain  per  pass  from  the 
200  pC  bunch  charge  can  tolerate  more  round-trip  losses  due  to  absorption  and  out- 
coupling. 

Extraction  for  a  Partially  Absorbant  Surface  with  a  Hole 


Figure  14.  Simulations  results  for  200  pC  bunch  charge  with  Q= 50 
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Figure  15  shows  the  a  comparison  between  the  optimal  hole  radius  of  0.05 
cm  for  100  pC  case  and  the  200  pC  optimal  hole  radius  of  0.15  cm  both  with  a  g~50. 
Although  the  dimensionless  amplitudes  are  not  very  different,  the  mode  shape  on  the  out- 
coupling  mirror  shows  that  this  set  of  parameters  will  include  higher-order  modes  as  the 
hole  size  increases. 


Figure  15.  Simulation  output  for  hole  radius  of  (a)  0.05  cm  with  100  pC  and  (b)  0.15  cm  with 

200  pC. 


c.  Varying  the  clipping  gap  for  the  theoretical  simulations 

Again,  simulations  were  conducted  varying  the  size  of  the  absorbing  gap 
over  the  range  of  0.3  cm  to  3.0cm. 

Figure  16  shows  the  results  of  the  simulations.  The  optimal  gap  in  this 
case  is  1.2  cm  with  an  extraction  of  3.0%.  After  the  peak  extraction  was  reached,  the  plot 
leveled  off  much  like  it  did  for  the  lOOpC  set. 
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Gap  Varitions  for  Bunch  Charge  of  200pC 


Figure  16.  Simulation  results  for  Gap  variations. 
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V.  MESH  OUT-COUPLING 


Mesh  out-coupling  is  another  out-coupling  technique  that  does  not  use  a 
conventional  partially  transparent  mirror.  Instead,  the  light  in  this  technique  will  out- 
couple  through  a  metallic  mesh.  There  are  two  types  of  meshes  that  can  be  considered 
for  this  application.  The  first  type  is  the  typical  metallic  mesh  that  uses  thin  strips  of 
metal  with  no  substrate  and  is  referred  to  as  an  inductive  mesh  [11].  The  second  type 
uses  a  thin  metallic  pattern  deposited  on  a  substrate  and  is  known  as  a  capacitive  mesh 
[1 1],  Both  can  have  either  circular  or  square  symmetry  and  are  illustrated  in  Figure  17. 
The  mesh  with  square  symmetry  can  be  parameterized  by  three  variables:  the  mesh 
period,  g,  the  grid  spacing,  a,  and  the  mesh  thickness,  t. 


Figure  17.  Illustrations  of  freestanding  (a)  inductive  metallic  mesh  and 
(b)  capacitive  mesh  with  square  symmetry. 

From  [1 1]  t  is  the  thickness  of  the  mesh,  g  is  the  mesh  period, 
and  a  is  the  grid  spacing. 

Mesh  out-coupling  has  been  proven  to  work  in  the  far  infrared  region  of  the 
spectrum  [12].  While  there  are  several  ways  to  implement  this  technique  [13],  two  will  be 
discussed  in  this  thesis:  single  mesh  out-coupling  and  double  mesh  out-coupling.  Both 
techniques  will  use  an  inductive  mesh  as  the  out-coupler  in  the  thin  mesh  regime  where 
t  «a ,  with  the  assumption  that  there  are  no  losses  due  to  absorption  in  the  mesh. 

A.  OPTICAL  PROPERTIES  OF  A  SINGLE  MESH 

The  optical  properties  of  a  single  mesh  can  be  derived  and  analyzed  to  show  that  a 
single  mesh  can  be  utilized  as  a  reflective  surface.  Previous  study  of  this  interaction  has 
provided  results  that  were  used  in  creating  parameters  for  this  thesis  [12].  The  analysis 


27 


from  this  study  showed  that  a  single  mesh  of  square  symmetry  has  optical  properties  that 
are  independent  of  polarization  because  the  mesh  can  be  considered  as  the  superposition 
of  two  crossed  strip  gratings  if  the  ratio  a/g<  0.2is  satisfied  [12].  The  same  paper 
argued  that  if  g/  A  <1,  i.e.  when  the  wavelength  is  longer  than  the  mesh  period, 
diffraction  from  the  mesh  edges  is  negligible  [12]. 


Many  of  the  coefficients  needed  to  describe  the  optical  properties  of  the  mesh  can 
be  calculated  based  on  an  equivalent  electric  circuit  that  yields: 

r  =  .  7/2 


1  +  7/2 
t=  i+r 


(V.l) 
(V.2) 

where  T  is  the  complex  amplitude  reflection  coefficient  of  the  mesh,  rt  is  the  complex 
transmission  coefficient  of  the  mesh,  Y  -  -i{2/ Z^Aj  g'j  is  the  admittance  of  the 
equivalent  circuit,  and  Z0  =  21n(csc(;ra/g)is  the  impedance  [12].  This  assumes  that 
there  are  no  losses  due  to  absorption  in  the  mesh  material. 


The  theoretical  power  reflectance  is  the  magnitude  of  the  complex  amplitude 
reflection  coefficient  of  the  mesh  squared,  R=\T\  .  This  can  be  used  to  show  how  the 
single  mesh  reacts  to  light  that  is  normally  incident  onto  it.  In  Figure  18,  the  power 
reflectance  was  plotted  as  a  function  of  g/A  and  parameterized  over  0.04  <a/g<  0.20  . 
The  figure  shows  that  when  X  is  much  larger  than  g,  the  single  mesh  tends  to  reflect  most 
of  the  power.  From  this,  it  can  be  seen  that  the  single  mesh  can  be  used  in  place  of  a 
mirror  in  the  optical  cavity  of  a  long  wavelength  FEL,  with  the  mesh  parameters  chosen 
to  produce  a  desired  quality  factor  for  the  cavity  for  a  fixed  wavelength. 
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Power  Reflectance  of  a  Single  Mesh  Out-Coupler 


Figure  18.  Power  reflectance  of  a  single  metallic  mesh  vs.  g/A  where  each  curve  is 

parameterized  by  a/g.  After  [12] 


B.  DOUBLE  MESH  OUT-COUPLING 

While  a  single  mesh  can  be  used  as  an  out-coupler,  a  more  flexible  design  uses  a 
double  mesh.  Instead  of  using  one  freestanding  metallic  mesh,  two  freestanding  metallic 
meshes  are  used  in  parallel  and  separated  by  some  distance  typically  on  the  order  of  an 
optical  wavelength,  A.  Then,  the  total  transmission  and  reflection  from  both  meshes 
depends  not  only  on  the  properties  of  each  individual  mesh,  but  also  by  the  separation 
between  them,  much  like  a  Fabry-Perot  Interferometer  (which  will  be  discussed  in  the 
following  section  of  this  chapter).  Whereas  the  optical  properties  of  the  single  mesh  were 
strongly  dependent  on  A,  the  extra  degree  of  freedom  afforded  by  a  double  mesh  allows 
one  to  maintain  the  constant  set  of  optical  properties  across  a  range  of  wavelengths  by 
adjusting  the  separation  between  the  meshes  [14]. 
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1.  Fabry-Perot  interferometer 

A  Fabry-Perot  interferometer  (FPI)  consists  of  two  partially  transparent,  parallel 
mirrors  separated  by  a  medium  or  vacuum.  The  FPI  is  also  used  in  a  variety  of 
applications  including  but  not  limited  to  optical  resonators  and  spectrometers. 

Figure  19  illustrates  how  a  FPI  transmits  and  reflects  light,  where  s  is  the  distance 
separating  the  two  parallel  reflecting  surfaces,  6  is  the  internal  angle  of  incidence,  and  E0 
is  the  initial  amplitude  of  the  incident  wave.  As  the  wave  travels  through  the  FPI,  each 
reflection  creates  a  partial  wave  that  is  either  transmitted  or  reflected.  The  subsequent 
internal  reflections  create,  in  principle,  an  infinite  number  of  partial  waves  [15], 


Figure  19.  Interaction  between  an  incident  light  beam  and  a  FPI. 


Each  partial  wave  acquires  a  phase  shift  per  pass  through  the  interferometer 
according  to 


Anns  cos  6 

A 


(V.3) 


where  n  is  the  index  of  refraction  of  the  medium  between  the  two  surfaces,  and  A  is  the 
wavelength  of  the  incident  wave,  Eq,i n  vacuum.  The  partial  waves  that  are  reflected  (Eo, 
E4,  Ee, . . .)  are  then  calculated  as 


:  rE{)  ,  E4  =  t2rE0e,s 


E6  =  t2r%e2iS 
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where  r  and  t  are  the  amplitude  reflection  coefficient  and  the  amplitude  transmission 
coefficient  respectively  and  are  assumed  to  be  equal  for  all  reflections  and  transmissions. 
Similarly,  the  partial  waves  that  are  transmitted,  {E\,  £3,  E5, ...),  are  also  calculated  as 

Ey  =  t2E0 ,  E3  =  t2r2E0eiS ,  E5  =  t2r4E0e2i\ 

The  sum  of  these  waves  constitute  the  amplitude  of  the  transmitted  wave  and  the 
reflected  wave, 


Et  +E3  +E5  +...+Elcr_x 

(V.4) 

Er  =E2+E4+E6+...+E2a 

(V.5) 

where  ET  and  ER  are  the  total  transmitted  amplitude  and  the  total  reflected  amplitude, 

respectively,  and  a  approaches  infinity.  When  the  values  of  the  partial  waves  are 

substituted  into  Equation  V.4  and  Equation  V.5,  they  become 

ET  =  t2(l  +  r2eiS  +  r4e2iS  +  ..)E0, 

(V.6) 

Er  =  [r  - +  t2reiS(  1  +  r2eiS  +  r4e2iS  +  ...)]^0 . 

(V.7) 

It  can  now  be  seen  that  the  total  complex  amplitudes  both  contain  an  infinite 
geometric  series  [15].  As  such,  the  amplitudes  can  be  rewritten  as 


T 

Et  "  \^E° 


(V.8) 
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1  -Re1 
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J-'n  5 


(V.9) 


where  R  =  r2  is  the  fraction  of  intensity  reflected  at  each  interface  and  T  =  t2  is  the 
fraction  of  intensity  transmitted  at  each  interface. 


If  the  intensity  of  initial  wave  is  defined  as  I0  -  E*0E0 ,  the  ratio  of  transmitted 
intensity  to  the  initial  intensity  can  be  described  as 

I_t  KEt  l1"*2) 

70  E*0E0  (l-i^  +  dtfsin2^/!)' 
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(V.10) 


Similarly,  the  ratio  of  reflected  intensity  to  the  initial  intensity  can  be  described  as 


T  F*  F 
1r  _  FrFr  _ 


4Rsin2  (8/ 2) 


(V.ll) 


(l-R2)  +  4Rsm2(S/2) 


An  assumption  of  Equations  V.8  through  V.ll  is  r2  +t2  =1,  i.e.  there  is  no  absorption 
occurring  in  the  mesh  material  [15]. 

Figure  20  illustrates  how  the  reflected  energy  depends  on  the  phase  shift  for  a 
given  R.  It  can  be  seen  that,  even  at  large  values  of  R,  there  are  phase  shifts  where  a  large 
amount  of  the  intensity  is  not  reflected.  From  equation  V.l  1,  when  8  =  2m7i,  where  in  is 
an  integer,  the  reflected  intensity  is  equal  to  zero.  As  such,  these  phase  shifts  allow  for 
complete  transmission,  assuming  that  there  is  no  absorption. 


Ratio  of  Reflected  Energy  to  Incident  Energy 


10 


5 


Figure  20.  Reflectance  of  a  single  metallic  mesh.  After  [15] 


Figure  21  displays  the  transmission  ratio  as  a  function  of  the  phase  shift.  As  seen, 
transmission  maximums  occur  for  the  same  phase  value  as  reflection  minimums. 
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Ratio  of  Transmitted  Energy  to  Incident  Energy 


Figure  21.  Transmittance  of  a  single  metallic  mesh.  After  [15] 

When  using  two  meshes  to  form  an  FPI,  the  only  change  to  the  FPI  equations  is 
the  reflectance,  R,  is  equated  to  the  square  of  T  from  Equation  V.l. 

2.  Spherical  Double  Meshes 

The  NPS  FEL  3D  code  assumes  that  the  optical  cavity  is  centered  about  the 
undulator  and  each  mirror  is  assigned  a  radius  of  curvature  based  on  the  Rayleigh  length 
of  the  cavity.  This  requirement  motivated  the  consideration  of  spherical  double  meshes 
to  ensure  the  feasibility  of  its  use  in  the  cavity  [16]. 

Suppose  a  light  beam  has  a  normal  incidence,  6  =  0 ,  on  the  interferometer  shown 
in  Figure  19.  In  that  case,  the  path  length  the  light  would  travel  through  the 
interferometer  is  2s  and  is  the  same  for  any  beam  normally  incident  along  the 
interferometer.  Therefore,  all  normally  incident  beams  exit  the  interferometer  with  the 
same  phase  change.  This  ensures  that  the  transmitted  and  reflected  beams  have  the  same 
phase.  If  a  set  of  spherical  surfaces  are  considered  for  the  interferometer,  it  must  be 
proven  that  the  path  length  traveled  is  approximately  uniform  for  any  normally  incident 
beam  at  any  point  on  the  surface  of  the  interferometer  in  order  to  preserve  the  uniform 
phase  change  of  the  reflected  and  transmitted  beams. 
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Consider  two  identical  spherical  meshes  that  are  separated  by  a  distance  s,  as  in 
Figure  22. 


Figure  22.  Beam  interaction  with  spherical  double  mesh.  Illustration  is  not  to  scale  as  z  ~ 

meters,  x  ~  centimeters,  and  s  ~  X. 

It  must  be  proven  that  there  is  no  appreciable  change  in  the  phase  of  the  beam  as 
it  travels  through  both  meshes  as  a  function  of  the  divergence  angle,  (p\  namely, 

Ac)  =  S (cp)  -  8 (0) «  2n  (V.12) 

where  AS  is  the  change  in  the  phase.  From  Equation  V.3,  it  is  shown  that, 
S ( (p )  =  (4 ns ')  j A  =  (4tt,s)/(/1  cos  (p)  and  S (0)  =  4 tts/ X .  When  substituted  into  Equation 
V.12,  the  result  is 

AS  =  [s/cos^-s]. 

Figure  22  shows  that  (p  =  tan  1  (x/z) « tan ‘(10  2)  « 10  2 rad.  From  this,  the  small  angle 

approximation  can  be  invoked,  cos  cp  «  cos(10“2) « 1 .  Figure  22  also  shows  that  s  ~  X. 
The  result  is  that  AS  «:  2n . 

3.  Equating  Quality  factor,  Q,  to  Mesh  Spacing 

In  order  to  simulate  the  out-coupling  interaction  of  the  double  mesh  out-coupler, 
the  metal  of  the  mesh  was  assumed  to  have  negligible  absorption  for  THz  radiation.  With 
this  assumption,  the  only  losses  are  transmission  through  the  mesh.  As  such,  the  quality 
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factor,  Q,  that  is  defined  to  be  inversely  proportional  to  the  loss  of  the  cavity,  is  also 
inversely  proportional  to  the  energy  transmitted  through  the  double  mesh,  Q^T~l .  If 
the  double  mesh  out-coupler  is  analyzed  as  an  equivalent  circuit  with  normalized 
admittances  of  each  mesh  being  Y\  and  Y2  [12],  the  power  transmittance  of  the  double 
mesh  out-coupler  can  also  be  calculated  in  the  following  fashion: 

rp  I  |2 

*  -|rzw|  ’ 

where  zDM  =  (cos(&s)-ry0  sin(fo,))(l  +  r/)  is  the  complex  amplitude  transmission 
coefficient  of  the  double  mesh  out-coupler  [12],  k  =  lx/ A  is  the  wave  number,  s  is  the 
spacing  between  the  meshes,  y0  =  (\  +  Y2  +  ztan(As))/(l  +  z(l  +  I^)tan(As))  is  the 

admittance  of  the  double  mesh  out-coupler,  T7  ={\-yI}l(\+yI)  is  the  amplitude 
reflectance,  and  y2  =  y^  +  Y^is,  the  input  admittance^  2],  From  the  previously  listed 
equations,  it  follows  that  a  relationship  between  Q  and  5  can  be  found,  as  shown  in  Figure 
23  (a).  The  periodic  nature  of  the  relationship  is  the  direct  result  of  the  periodic  nature  of 
the  interference  effects  between  the  meshes,  which  can  be  seen  in  Figure  23  (b)  where  the 
reflectance,  R=\-T,  is  plotted  versus  equivalent  mesh  spacing. 
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Figure  23.  (a)  Quality  factor  as  a  function  of  mesh  spacing. 

(b)  Reflectance  as  a  function  of  mesh  spacing. 
Oscillations  occur  with  a  period  of  rik! 2. 
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The  two  plots  in  Figure  23  also  specify  a  range  of  values  for  the  mesh  spacing. 
From  graphical  analysis,  5  is  predicted  to  have  a  usable  range  of  value  between  96pm  <  s 
<  103pm.  Figure  23  (b)  shows  that  if  5  is  less  than  96pm,  the  reflectance  approaches  zero 
and  that  would  not  allow  net  gain  in  the  optical  cavity.  If  5  is  greater  than  103  pm,  the 
reflectance  approaches  unity  at  which  point  no  light  is  being  out-coupled  from  the  cavity. 
Due  to  the  cyclic  nature  of  the  interferometer,  there  are  also  other  values  of  5  that  yield 
the  same  optical  properties.  In  this  thesis,  only  the  values  around  100pm  will  be 
considered. 

The  motivating  factor  for  finding  an  equivalent  mesh  spacing  for  a  given  Q,  flows 
from  the  prohibitive  nature  of  creating  a  computational  model  of  the  mesh  structure.  The 
difficulties  lie  in  that  the  model  would  need  to  be  on  the  scale  of  centimeters  with 
openings  on  the  scale  of  microns.  And  each  opening  would  need  to  have  enough  grid 
points  over  it  to  satisfactorily  represent  the  results.  This  would  require  large  amounts  of 
computational  time,  on  the  order  of  weeks  per  simulation  run.  In  the  approximations 
used  here,  the  macroscopic  optical  properties  of  the  double  mesh  can  be  used  with  a 
degree  of  confidence  in  the  results. 

C.  SIMULATION  RESULTS 

With  the  relationship  between  Q  and  5  established,  simulations  were  conducted 
using  the  NPS  3D  oscillator  FEL  code  and  the  parameters  listed  previously,  and  again 
categorized  into  two  groups  based  on  the  electron  bunch  charge. 

1.  lOOpC  Bunch  Charge  Simulations 

Simulations  were  conducted  varying  the  quality  factor  over  the  range  5<  Q  <100 
in  the  simulation  code,  to  represent  a  corresponding  mesh  spacing,  s,  as  discussed  in  a 
previous  section.  The  mesh  parameters  used  in  the  conversion  of  Q  to  s  were  a/g=O.  I6 
and  g/A= 0.20  which  corresponded  to  a  single  mesh  power  reflectivity  of  95%.  These 
parameters  satisfy  the  condition  for  long  wavelengths. 

Figure  24  shows  that  as  the  spacing  is  increased,  the  extraction  increases  to  a  peak 

and  appears  to  begin  to  decrease.  If  the  reflectance  is  too  low,  there  is  too  little  net  gain 

in  the  cavity  to  exceed  the  losses.  If  the  reflectance  is  too  high,  the  fields  inside  the 
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cavity  build  up  to  such  a  large  value  that  the  electrons  overbunch,  reducing  the  extraction, 
as  discussed  in  the  previous  chapter..  The  extraction  plot  also  shows  that  there  is  a  peak 
extraction  value  of  2.88%  when  R& 98% which  equates  to  Q=75,  and  an  equivalent 
mesh  spacing  of  s  ~  101/rni . 


Figure  24.  The  extraction  from  a  double  mesh  out-coupler  with  bunch  charge 
of  100  pC.  The  Equivalent  Mesh  spacing  directly  corresponds 
to  the  range  of  Q :  5<Q<95. 

Figure  25  displays  the  output  diagrams  for  the  peak  extraction.  It  should  be  noted 
that  the  mode  is  close  to  the  fundamental  Gaussian  as  seen  in  Figure  25  (b).  This  is 
preferable  because  the  energy  produced  can  then  be  focused  into  a  smaller  spot  in  an 
application. 
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Figure  25.  (a)  Simulation  output  of  the  evolution  of  the  optical  field 

inside  the  undulator.  (b)  Optical  mode  shape 
at  the  out-coupling  mesh  surface. 


2.  200  pC  Bunch  Charge  Simulations 

Simulations  were  also  conducted  for  a  bunch  charge  of  200  pC  by  varying  the 
quality  factor  over  the  range  of  Q  values  previously  stated.  Figure  26  illustrates  how  the 
simulated  FEL  performed  with  the  increase  in  bunch  charge.  The  plot  shows  that  the 
peak  extraction,  2.9%  was  approximately  the  same  as  the  100  pC  case.  The  optimal 
reflectivity  was  R= 0.97,  which  corresponds  to  a  mesh  spacing  of  5  «  99jum  and  Q= 35. 
However,  the  simulated  FEL  did  perform  better  over  a  broader  range  of  mesh  spacing. 
This  can  be  attributed  to  the  fact  that  the  increase  in  the  bunch  charge  directly  correlates 
to  an  increase  in  net  gain  per  pass. 
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Figure  26.  The  extraction  from  a  double  mesh  out-coupler  with  bunch  charge 
of  200  pC.  The  green  dotted  line  represents  the  results 
from  the  100  pC  case.  The  Equivalent  Mesh  spacing  directly 
corresponds  to  the  following  range  of  Q :  5<Q<95. 

Figure  23  displays  the  output  diagrams  for  the  peak  extraction.  It  should  be  noted 
that  the  mode  shape  in  Figure  27  (b)  is  once  again  close  to  the  fundamental  Gaussian. 


Figure  27.  (a)  Simulation  output  of  the  evolution  of  the  optical  field  inside  the 

undulator  field  amplitude  in  the  optical  cavity. 

(b)  Optical  mode  shape  at  the  out-coupling  mesh  surface. 
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3. 


Conclusion 


The  double  mesh  out-coupler  simulations  showed  the  FEL  had  extraction  levels 
comparable  to  hole  out-coupling  and  traditional  partially  transparent  mirrors,  as  would  be 
expected.  Unlike  hole  out-couplers,  mesh  out-couplers  will  not  disturb  the  transverse 
mode  shape  as  long  as  X  »  g.  For  larger  values  of  mesh  spacing,  the  100  pC  bunch 
charge  may  have  a  larger  range  of  relatively  high  extraction,  whereas  the  200  pC  bunch 
charge  has  better  performance  even  at  small  spacing  values. 
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VI.  FUTURE  WORK 


Several  THz  FEL  out-coupling  methods  were  studied  in  this  thesis,  including  hole 
and  double  mesh  out-coupling.  Another  method  uses  zone  plate  diffractive  optical 
elements.  Zones  plates  are  surfaces  that  use  diffraction  to  focus  incident  light.  As  seen 
in  Figure  28,  the  zone  plate  is  comprised  of  concentric  rings,  or  zones,  that  alternate 
between  being  transparent  and  opaque.  The  transparent  zones  are  spaced  to  promote 
constructive  interference  of  the  light  at  a  specified  focal  point.  To  achieve  this,  the  zone 

radii  must  be  determined  by  rn  =  JnA /  +  (n2/ 12  JJ4  ,  where  n= 1,  2,  3,  ...,  is  the  integer 

number  of  the  zone,  /  is  the  distance  to  the  focus  point,  and  X  is  the  wavelength  of  light 
being  focused  [17], 


Figure  28.  Zone  plate  with  5  zones.  Blue  areas  are  opaque  and 

red  areas  are  transparent. 


Figure  29  shows  the  diffraction  pattern  in  the  plane  of  a  zone  plate  at  the  focus 
resulting  from  five  zones.  Constructive  interference  produces  a  small,  central  peak  at  the 
focus,  with  a  faint  series  of  concentric  rings  due  to  higher  diffraction  orders. 
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Figure  29.  (a)  2D  output  intensity  of  light  out-coupled  using  a  zone  plate  at  focal  point. 

(b)  3D  output  intensity  of  light  out-coupled  through  a  zone  plate 


This  unique  method  for  out-coupling  was  previously  considered  to  analyze 
whether  it  was  possible  to  use  a  zone  plate  in  a  resonator  [18].  The  analysis  showed  that 
under  certain  conditions,  namely,  when  the  focal  length  is  less  than  the  resonator  length, 
a  zoned  out-coupler  would  be  feasible  [18],  [19]. 

The  NPS  3D  FEL  code  is  currently  being  modified  to  simulate  mirrors  with  zone 
structures.  When  this  modification  is  complete,  simulations  should  be  conducted  in  order 
to  examine  the  viability  of  this  method  for  use  in  the  THz  regime. 
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VII.  CONCLUSION 


There  are  many  current  and  potential  applications  of  intense,  coherent  THz  light 
sources,  including  basic  physics  experiments,  military  and  homeland  security 
applications.  The  FEL  is  already  a  powerful  source  of  coherent  light  in  many  parts  of  the 
electromagnetic  spectrum,  including  THz.  From  this  motivation,  this  thesis  explored 
multiple  out-coupling  techniques  for  an  FEL  operating  ~  100pm. 

Through  simulation  using  the  NPS  3D  FEL  code,  hole  out-coupling  was  shown  to 
be  a  viable  method  for  out-coupling  in  the  THz  regime.  It  was  seen  that  the  hole  had  an 
effect  on  the  optical  mode  in  the  cavity,  and  that  there  was  an  optimal  hole  size  for 
maximum  extraction.  This  optimal  hole  size  was  shown  to  be  affected  by  changes  in  the 
bunch  charge  and  the  amount  of  absorption  of  the  mirror  material.  The  hole  was  also 
shown  to  affect  the  transverse  mode  shape  of  the  beam;  when  the  hole  radius  was 
increased,  alterations  in  the  mode  shape  increased. 

As  suggested  from  a  previous  study,  beam  clipping  at  the  ends  of  the  undulator 
was  shown  to  improve  mode  quality  when  using  hole  out-coupling.  An  optimal  amount 
of  clipping  was  determined  for  the  optimal  hole  radii  studied  in  this  thesis. 

Double  mesh  out-coupling  was  also  shown  to  be  potentially  a  good  method  for 
THz  out-coupling.  The  reflectance  of  the  double  mesh  out-coupler  was  shown  to 
increase  as  the  spacing  between  the  meshes  increased.  An  optimal  reflectance, 
corresponding  to  an  optimal  mesh  spacing,  was  determined  for  different  values  of  bunch 
charge.  It  was  also  shown  that  mesh  out-coupling  has  very  little  effect  on  the  optical 
mode  shape.  A  double  mesh  has  the  additional  advantage  that  it  can  be  adjusted  to 
perform  well  over  a  broad  range  of  wavelengths,  to  take  advantage  of  the  FEL’s  inherent 
tunability. 
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